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ABSTRACT: Employing hamiltonians denned by two-body embedded Gaussian 
orthogonal ensemble of random matrices(EG0E(2)) plus a mean-field producing one- 
body part, strength functions (for states defined by the one-body part) are constructed 
for various values of the strength of the chaos generating two-body part. Numerical 
calculations for six and seven fermion systems clearly demonstrate Breit-Wigner to 
Gaussian transition, in the chaotic domain, in strength functions as found earlier in 
nuclear shell model and Lipkin-Meshkov-Glick model calculations. 



PACS number(s): 21.10.Pc, 24.60.Lz, 21.60.Cs 



Strength functions (also called local density of states (LDOS) in literature) of 
simple modes are basic ingredients of many particle systems such as atomic nuclei 
|Ij]. In the last few years, with developments in quantum chaos 0, nature of strength 
functions for finite isolated interacting many particle systems are being investigated 
using a variety of models: (i) Zelevinsky et al using nuclear shell model ||; (ii) Wang et 
al using the three-orbital Lipkin-Meshkov-Glick model (iii) Borgonovi et al using 
a symmetrized coupled two-rotor model ||; (iv) Benet et al using a chaotic model of 
two coupled quartic oscillators ||. Results from the shell model and Lipkin-Meshkov- 
Glick model studies clearly showed that strength functions exhibit, for interacting 
particle systems, Breit-Wigner to Gaussian transition. Moreover it is seen that the 
transition takes place much after the onset of chaos in energy levels (i.e. much after 
level fluctuations start following GOE). In order to establish that this transition is 
generic to interacting particle systems, we carried out one plus two-body embedded 
Gaussian orthogonal ensemble (EGOE(l+2)) calculations and the results are reported 
in this short communication (see also |7L M). 

Given a compound state <f>k, the probability of its decay into stationary states ipE 
(generated by H) is given by \{4>k I ^e)^- Then the strength function F k (E) is, 

l&>=£C*IVto> => F k (E) = Y: \cf\ 2 5(E - E>) = (5(H - E)) k (1) 
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The standard form ||, normally employed in many applications in nuclear physics, 
for strength functions is the Breit-Wigner (BW) form characterized by a spreading 
width Tfc, 
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F k ..sw{E) = - {E _ m k 2 + n/A (2) 

where E~ k = (k \ H \ k) = /f^ F k (E)E dE. Starting with H = h(l) + XV (2) (a one 
plus two-body interaction as in the shell model) with h(l) defining the <pk states, it 
is easily seen that the assumptions that give Fk-.Bw(F) will break down when the 
mixing is strong. Thus it is expected that the form of F k {E) will be different from 



BW for large A. In order to find the form of strength functions in large A limit, it 
is useful to think of (f>k as a compound state generated by the action of a transition 
operator O on a state ipE k (for example ground state), 
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((OW5(H - E k ))) 

The numerator in the expression for F^{E) in (3) is nothing but the bivariate strength 
density generated by the operator and then it is easy to recognize that Fk(E) 
is the conditional density of this bivariate strength density It is known || that 



EGOE(/c), the embedded GOE of /c-body interactions (see Ref. (K| for the definition 
of EGOE(k)), in general gives bivariate Gaussian strength densities. Therefore, for 
EGOE(2), which models a generic two-body hamiltonian, the strength functions take 
Gaussian form (conditional density of a bivariate Gaussian is a Gaussian). However, 
for generic interacting particle systems, it is more appropriate to consider one plus 
two-body EGOE(l+2) hamiltonians {H} = h(l) + A{V(2)} where h(l) is the one- 
body mean-field producing part and {V(2)} is the chaos generating two-body part; 
note that {V(2)} is GOE in two-particle space with unit matrix elements variance 
and A is the interaction strength. Starting with EGOE(l+2) it is to be expected that 
for sufficiently large values of A, EGOE(2) description should be valid and therefore 
(3) gives the shape of the strength function to be Gaussian for large A; this result 
is indeed seen in numerical calculations. Now the important questions are: (i) how 
Fk(E) changes as A is varied and (ii) from which value of A strength functions take 
Gaussian form. 

For further understanding of the nature of Ff-(E), we performed EGOE(l+2) 
calculations in 924 dimensional space generated by six fermions (m = 6) in twelve 
single particle states (N = 12) and similarly in the 3432 dimensional iV = 14 and 
m = 7 space. The single particle energies employed are e, = i = 1, 2, . . . , iV as 



in [11]. For various values of A in h(l)+X{V(2)}, strength functions are constructed by 
choosing the (f>k states to be the mean-field h(l) states defined by the distribution of m 
particles in the N single particle states; their energies Sk are Sk = (k\h(l) + X V(2)\k). 
In the calculations E and Sk are zero centered for each member and scaled by the 
spectrum (E'a) width a; = (Sk — e)jo and E = (E — e)/a. For each member \C^\ 2 
are summed over the basis states \k) in the energy window E}~±A and then ensemble 
averaged (E) vs E is constructed as a histogram; the value of A is chosen to be 
0.025 for A < 0.1 and beyond this A = 0.1. Results for E\~ — are shown for A = 
0.05, 0.1, 0.15, 0.2, 0.3 and 0.5 in Fig. 1; in the plots J F £k (E) dE = 1. From the 
figures it is clearly seen that there is BW to Gaussian transition in Fk(E). A measure 
for this transition is, 

RW = ^7 T2 (4) 

i 

In (4) Ei are defined by the center of each bin in the histogram representing Fj, x \E) 
and ED corresponds to Gaussian with Edgeworth corrections; the ED incorporates 
12] skewness (71) and excess (72) corrections. Strength function F k (E) is BW for 
R = and Gaussian (or ED) for R = 1. The value of the interpolating parameter 
A = Xp k for onset of the transition from BW to Gaussian is taken to be i?(A^ fc ) = 0.7. 
For A = 0.05, 0.1, 0.15, 0.2, 0.3, 0.5 the values of R are 0.13, 0.06, 0.21, 0.7, 0.9, 
1.03 for the six particle example. For the seven particle example the corresponding 
numbers are 0.15, 0.02, 0.45, 0.75, 0.93, 1.03 respectively. Thus, with the measure R 
in (4), Xp k « 0.2 for both the 6 and 7 particle cases. 

In Fig. 2 shown are EGOE(l+2) results for various values of A for the nearest 
neighbour spacing dustribution and A 3 (L), < L < 40 statistic. The \ Fk deduced 
from the results in Fig. 1 should be compared with A = A c ~ 0.08 and 0.06, derived 
via Fig. 2, for six and seven particle examples repectively for onset of chaos in level 



fluctuations. The A c numbers are consistant with (Nm 2 ) 1 variation as expected 
from the theory given in | 13| . It should be noted that in the (2s, l c Q m=12 >- rrT=0+0 



shell model example considered in ||, A c ~ 0.3 and Xp k ~ 0.6. Thus the BW form 
for Fk{E), which begins some what before A approches A c (for A << A c the strength 
functions are delta functions with perturbative corrections), extends much into the 
chaotic domain (defined by A > A c ) and the transition to Gaussian shape takes place 
in the second layer defined by A^ fc (A^ >> A c ). As the A^ fe in our two examples did 
not show much variation with (N, m), it appears that \F k will have weak dependence 
on m unlike A c . Perhaps for sufficiently large m, the transition to Gaussian form takes 
place in the thermalization regime discussed by Flambaum and Izrailev |7|] (see also 
HU). It should be remarked that EGOE(2) gives Gaussian form for state densities 



for sufficiently large m and its extension to partial state densities (which are sums of 



strength functions) is indeed the basis of statistical nuclear spectrocopy || [15 



In conclusion, the results of the present EGOE(l+2) study (Figs. 1,2) and the 
earlier shell model || and Lipkin-Meshkov-Glick model j|] analysis establish firmly 
the Breit-Wigner to Gaussian transition for strength functions, in the chaotic domain, 
in interacting particle systems such as atomic nuclei. 
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FIGURE CAPTIONS 



Fig. 1 Strength functions for EG0E(l+2) for various values of the interaction 
strength A: (i) for a system of 6 fermions in 12 single particle states with 25 mem- 
bers; (ii) for a system of 7 fermions in 14 single particle states (due to computational 
constraints, here only one member is considered just as in |TT|). In the figure, the 
histograms are EGOE(l+2) results and continuous curves are BW fit. For the 6 
fermions case, the dotted curves are Gaussian for A < 0.15 and Edgeworth corrected 
Gaussians (ED) for A > 0.15. Similarly for the 7 fermions case, the dotted curves are 
Gaussian for A < 0.1 and Edgeworth corrected Gaussians (ED) for A > 0.1. See text 
for further details. 

Fig. 2 Nearest neighbour spacing distribution P(S) vs S and Dyson-Mehta As(L) 
statistic for < L < 40 for various values of the interaction strength A in EGOE(l+2). 
For P{S), histograms are EGOE(l+2) results, dashed curves are Poisson and contin- 
uous curves are Wigner distribution. For A 3 (L), filled circles are EGOE(l+2) results, 
dashed curves are Poisson and continuous curves are for GOE. Results are shown for 
the six and seven particle examples considered in Fig. 1. 
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